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We study massive spin-1/2 fermion production in nonsingular superstring cosmology, taking into 
account one-loop quantum corrections to a superstring effective action with dilaton and modulus 
fields. While no creation occurs in the massless limit, massive fermions can be produced by the 
existence of a time-dependent frequency. Due to the increase of the Hubble expansion rate during 
the modulus-driven phase, the occupation of number of fermions continues to grow until the point 
of the graceful exit, after which fermion creation ceases with the decrease of the Hubble rate. 
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I. INTRODUCTION 

With the development of superstring theory, string- 
inspired cosmological models jl| have received much at- 
tention to describe the evolution of the very early stage 
of the Universe. Most of such scenarios are based on the 
low-energy effective field theory, which is expected to be 
valid at the Planck scale. While a full theory is not yet 
established, it is important to test the viability of string 
theories by extracting various cosmological implications 
from them. 

Among string-motivated cosmological models pro- 
posed so far, the pre-big-bang (PBB) scenario has been 
most widely studied. If one assumes that the Universe 
has a T-duality, there exist two disconnected branches. 
One of which (t < 0) corresponds to the stage of super- 
inflation driven by the kinetic term of the dilaton field, 
and another (t > 0) is the Friedmann branch where the 
Universe exhibits standard decelerating expansion. The 
PBB scenario basically requires the existence of nonsin- 
gular solutions which interpolates between two discon- 
nected branches ||. In the tree- level superstring action, 
however, one has no-go results that singularity can not 
be avoided [|J|. 

In order to overcome such singularity problems, Anto- 
niadis, Rizos, and Tamvakis involved one-loop quan- 
tum corrections to the string effective action with dila- 
ton and modulus fields, and found some nonsingular 
solutions. Since the success of singularity avoidance 
is mainly determined by the motion of the modulus 
field, the allowed ranges of parameters have been an- 
alyzed in the absence of the dilaton field in the flat 
Friedmann-Robertson- Walker (FRW) background j?J and 
the anisotropic Bianchi type-I metric In the full 

dilaton-modulus system, several authors studied gener- 
ality of singularity avoidance in the closed FRW Q and 
Bianchi type-I and -IX |l0| spacetimes. It was found 
that nonsingular solutions generically exist except for the 
Bianchi-IX case. 

From observational point of view, analysis of the per- 
turbations predicted in the PBB model is an important 
issue in order to test realistic string theories. In this 
respect, many authors investigated quantum creation of 



scalar particles such as dilaton and axions 
production of gravitational waves [p"7|-|l9| , most of which 
exhibit different spectra compared to the standard cos- 
mology. Another interesting prediction in the PBB sce- 
nario is the generation of primordial magnetic fields due 
to the break of the conformal invariance (2Cjj2l|] . 

Recently, Brustein and Hadad [^2| studied fermion 
production in superstring cosmology in the presence of 
the dilaton coupling only. It was found that mass- 
less fermions are not created, since the equation of 
fermions reduces to that in the Minkowski spacetime in 
the massless limit. In this Letter, we investigate massive 
fermion production in more general dilaton-modulus sys- 
tem whereby singularity can be avoided. In fact we will 
show that massive fermions are nonadiabatically created 
due to the existence of a time-dependent mass term. 



II. NONSINGULAR SOLUTIONS 

Consider the following one-loop effective action of the 
heterotic superstring theory |6PJIo[] , 
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written in the Einstein frame. Here R, 4>, a, and H^x 
are the scalar curvature, the dilaton, the modulus, and 
the antisymmetric tensor field, respectively. In this work 
we set H^ux — and neglect the curvature terms higher 
than the second order. The Gauss-Bonnet term, Rqb, is 
defined as 
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In the presence of the last term in the action ( |2.l[ ) (i.c, 
one-loop quantum corrections), singularity problems in 
the tree-level action can be avoided ||. The coefficients, 
A and S, are determined by the inverse string tension a' 
and the four-dimensional trace anomaly of the N = 2 
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sector, respectively. While A is positive definite, S can be 
either positive or negative. 

The function, £(<t), is expressed as 
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Then the first derivative of £(cr) in terms of a is well 
approximated as £'(cr) = — (27r/3) sinhcr, which we use 
in our numerical analysis. 

It is also convenient to introduce a dimensionless func- 
tion, f((f>, a) = [e* - 6£{a)]/16 with 5 = S/X. We nor- 
malize time and spatial coordinates by the string length 
scale y/~X as x^ = x 11 /\f~X 1 and scalar fields as <fi — <f>VX, 
a = <t\/A. Hereafter we drop bars for simplicity. 

Adopting the flat FRW metric as the background 
spacctimc, with a = e p being the scale factor, the dy- 
namical equations for the metric and scalar fields yield 



+ 8pf)(p + p 2 ) + 4(1 + 8/)p 2 + cf 1 + 3a 2 
+ 3p0 - 2f^R 2 GB = 0, 
2 
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together with the constraint equation, 
Up 2 + 96p 3 / - ft - 3ct 2 = 0. 



(2.4) 
(2.5) 

(2.6) 
(2.7) 



Here an overdot denotes a derivative with respect to 
cosmic time, t, and the Gauss-Bonnet term is given as 
Rq B = 24p 2 (p + p 2 ). Nonsingular cosmological solutions 
have been found for negative values of 6 
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FIG. 1. The evolution of the Hubble rate is plotted when 
the singularity is avoided. We choose <5 = — 48/-7T, and set 
4> = 4> — 0, a — 0, p = at t = 0. a is determined by 
the constraint equation (2.7) as a = 0.2. Inset: (f> and a vs 
t. The evolution of the system is dominated by the a field 
around the graceful exit. 

In the absence of the modulus field, the scale factor 
in the Einstein frame evolves as a a It] 1 / 3 during the 
dilaton-driven phase. For negative i, this corresponds 



to the accelerated contraction, d < and a < 0. In 
the tree-level action, one needs to assume that the epoch 
of the accelerated evolution comes to an end at some 
time in order to make a smooth transition to another 
branch (t > 0). In the present scenario, however, taking 
into account one-loop corrections opens up the possibility 
of the graceful exit driven by the kinetic energy of the 
modulus field. 

We show one nonsingular solution in Fig. 1. Generally 
the Hubble rate, H = p, grows as H cx (— 1)~ 2 during 
the modulus-driven phase. In order to avoid singularity 
at t = 0, the velocity of (f> is required to be much smaller 
than a during the graceful exit. If \§\ is sufficiently small 
around t = 0, there always exist nonsingular solutions for 
4> < 0. This is because « for negative large values of 
(f>, which indicates that the singularity avoidance is prac- 
tically independent of <p. The allowed parameter regions 
with respect to 4> were precisely analyzed in Ref. [ pi , 
which we do not repeat it here. 

The sign of & does not change during the whole evolu- 
tion. For a > which corresponds to the case of Fig. 1, a 
rapidly moves relative to <p around the graceful exit. Af- 
ter a smooth transition at t = 0, the dilaton freezes with 
<j) ~ and the modulus evolves as a cx In t as t — > oo. 
The Hubble rate begins to decrease for t > 0, which 
asymptotically approaches the Friedmann-like Universe, 
a cx i 1 / 3 and H - l/(3i). 



III. MASSIVE FERMION PRODUCTION 

Let us consider the following action for the fermion 
field if) with bare mass m: 
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d 4 xy/=g f (iip^Tl^ - mtpt/;) , (3.1) 



where 7 M is the curved-space Dirac matrices, and = 
+ (l/4)7 Q/ 3o;^ /3 is the spin-1/2 covariant derivative, 
where wff is the spin connection. 7 Q denotes the Dirac 
matrices in Minkowski spacetime with "f a p = "f[ a J/3] ■ 
From the action (13 . If) we obtain the Dirac equation, 



0. 



(3.2) 



Since 7 =7°, 7* 
in the flat FRW background, Eq. (3.2) is simplified by 
introducing a new field, \ = a 3/>2 f 1 ^ 2 ^, as 



7 7a, and 7^ = 7^ + (3/2)ff 7 ° 



{il^dfj, - ma) \ = 0, 



(3.3) 



where 8q denotes the derivative with respect to conformal 
time, rj = j a _1 iit. We decompose the x field into Fourier 
modes as 
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[u,(fr, il)a s {k) + v s (k, vHi-k)] , (3.4) 
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where v s (k) = Cu^(—k) with C being a constant. 

Defining u s = [u+(r])tp s (k), su-(r))cp s (k)] T and v s — 
[sv+(r])ip s (k),V-(r])ip s (k)] with ip s (k) being eigenvec- 
tors of helicity operators, the Dirac equation (3.3) reads 

MM 



u'±(rj) = iku T (r]) =p imau±(rf), 
which reduces to the decoupled form: 
u± + [uil ± i(ma)'\ u± = 0, 



(3.5) 



(3.6) 



where uJf, = k 2 + (ma) 2 . Note that we imposed 
the normalization conditions, u\(k : r])v s (k : r]) = 0, 
ul(k,-q)u s (k,ri) = v\ (k, rf)v s (k, rf) = 5 rs , \u + \ 2 + \u-\ 2 = 
2. 

In order to diagonalize the Hamiltonian, we introduce 
new operators, a(k,r]) — ctk(r])a(k) + /3fc(?7)^(— k) and 
tt(k,r)) = -0l(ri)a(k) + a^(?7)6 t (-fc), where the Bo- 
golyubov coefficients satisfy 



"fc = — ™ — Pk, \pk\ 
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with 



E k = fcRe(u+u_) + ma (l - |u+| 2 * 
F k = (u 2 + — it_J + mau + U-. 



(3.7) 



(3.8) 
(3.9) 



Note that the canonical commutation relation leads to 
Wk\ 2 + |/3fe| 2 = 1, which restricts the occupation numbers 
of fermions as n k = \(3k\ 2 < 1- The initial conditions are 
chosen as u±(r)o) — (a^ T ma)/u>k, corresponding to 
Ek(rio) = Wfc, Fk(r)o) = 0, and nk(r]o) = 0. 



1" 1 

r 


1 


r 

r k = 10" 2 , 


m = 10 " 2 / k = 10 1 ,m = 13 

^/ A/ i 




, 






F ^ . ' 

F k = l,m = 10" 6 1 



1.0 
0.010 
0.00010 
10"' 
10' 8 
* 10 10 

io- 12 
io- 14 

10" 

io- 18 

10 .20 



FIG. 2. The evolution of the occupation number of 
fermions for three cases: k = lCP 2 ,m = 10~ 2 ; 
k = 10 _1 ,ra = 1; and k — l,m = 10 -6 , where k and m are 
normalized by the string length scale, Va. The S and initial 
values of background quantities are the same as in Fig. 1. The 
enhancement of fermions strongly depends on the relation of 
two parameters, k and m. 



In the limit of m - 
Minkowski spacetime, 

«4 



0, Eq. (3.6) reduces to that in 



k 2 



u± = 0. 



(3.10) 



The solution for this equation is expressed as u±(rj) 
where we used the initia l co n ditio ns, u±(rjo) 



1. Then we have nk = by Eqs. ( |3.7[ )-(3.9), which indi- 
cates that no creation occurs in the massless limit p2]| . 

When the mass of fermion is sufficiently small relative 
to the physical wave number (to -C k/a), the situation 
is similar to the massless case. In Fig. 2 we find that 
fermions are hardly excited for k — 1 and m = 10~ 6 , 
where k and to are normalized by the string length scale, 
Note that the acquired number of e-foldings during 
the modulus-driven phase is not large compared to the 
standard inflationary scenarios, e.g., if we normalize the 
scale factor as a = 1 for t = 0, a w 0.345 for t = —500 
at which the contribution of modulus begins to be im- 
portant relative to dilaton. This indicates that when the 
condition, to <C k/a, holds in the initial stage of the 
modulus-driven phase, it is typically valid even around 
the graceful exit. 
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FIG. 3. The occupation number of fermions at the end 
of the modulus-driven phase as a function of the fermion 
mass, m, for three different comoving momenta, k — 10 -4 , 
k = 1CP 2 , and k = 1. The 8 and initial values of background 
quantities are the same as in Fig. 1. 

When to and k/a are comparable during the modulus- 
driven phase, massive fermions are created due to the 
existence of the time-dependent mass term in Eq. 
In Fig. 2 the number density of fermions continues to 
grow until the point of the graceful exit for the cases of 
k = 10~ 2 ,to = 10 -2 ; and k — IO" 1 ,?™ = 1. The nona- 
diabatic condition where particles are sufficiently excited 
can be written as \tOk\ > w 2 , which yields 

[fc 2 + (TOa) 2l3/2 



H > 



(ma) 2 



(3.11) 



In a dust or radiation dominated Universe, the Hubble 
rate decreases as H cx 1/t, which works to violate the 
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nonadiabatic condition, (3.11). In fact, in preheating af- 
ter inflation, unless an inflaton decay to fermions is not 
taken into account, the mass term (ma) 2 does not lead 
to sufficient fermion production p3| , p5| -p7| . However, in 
the present model, the growth of the Hubble rate during 
the modulus-driven phase assists the nonadiabatic con- 
dition to hold, which results in nonpcrturbative particle 
creation solely by the time-dependent mass term. The 
growth of the occupation number ends after the smooth 
transition to the Friedmann-like Universe, since the Hub- 
ble rate begins to decrease (see Fig. 1). 

When to is much larger than k/a, fermion production 
is generally suppressed. Espe cially for to 3> H wher e the 
nonadiabatic condition (3.11) is not satisfied, Eq. ( |3.6| ) 
is approximately written as X± + m 2 X± ~ with X± = 
a 1 ' 2 u±. This indicates that fermions are hardly created 
in the supermassive limit, to — > oo. 

In Fig. 3 we plot the occupation number of fermions 
at t = as a function of mass for three different mo- 
menta. For each momentum, there exists a maximum 
rife for some value of to (= to*). Since to, is typically 
of the same order as the each corresponding momentum, 
the curves shift from left to right with increasing k. If 
the 1/vA is around the Planck scale, our results suggest 
that massive fermions heavier than the GUT scale can 
be copiously produced, which may play important roles 
for the leptogenesis scenarios EJ . 



to zero depending on two coupling constants, h\ and /12, 
this may further strengthen nonadiabatic amplification 
of fermions. We leave to future work about the precise 
investigation of this issue. 

Although we have restricted ourselves in spin-1/2 
fermions satisfying the Dirac equation, nonthermal pro- 
duction of gravitinos (spin-3/2 fermions) has recently be- 
come an issue of great importance pjjj] . Gravitinos have 
both helicity-3/2 and -1/2 states. While the helicity- 
3/2 mode reduces to the form of the Dirac equation, the 
helicity-1/2 mode behaves like the goldstino in global su- 
persymmetric limit. It was shown in Ref. ]22| that the 
latter mode also reduces to the same form as the former 
mode in the massless limit by assuming the power-law 
evolution of the scale factor in the standard PBB sce- 
nario, which means that creation of massless gravitinos is 
highly suppressed. However, the situation will change for 
massive gravitinos, in which case the existence of time- 
dependent mass terms may lead to the amplification of 
gravitinos during the modulus-driven phase. It is cer- 
tainly of interest to study gravitino production in real- 
istic nonsingular cosmology, since this provides a pow- 
erful mechanism to distinguish between different models 
of string theories and regions of the parameter space, 
together with the CMB constraint by scalar and tensor 
metric perturbations produced during the graceful exit 

M. 



IV. CONCLUSIONS AND DISCUSSIONS 
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We have investigated the production of massive spin- 
1/2 fermions in nonsingular superstring cosmology with 
dilaton and modulus fields. The existence of the modulus 
coupled to the Gauss-Bonnet curvature invariant leads 
to a smooth transition from the modulus-driven accel- 
erated expansion phase toward the Friedmann-like Uni- 
verse. Since the Hubble rate increases before the graceful 
exit, this makes it possible to produce massive fermions 
in terms of nonadiabatic change of their frequencies. In 
particular, fermions are most efficiently excited when the 
bare mass to is comparable to the physical momentum, 
k/a. In both massless and supermassive limits, fermion 
creation is strongly suppressed. 

The occupation number of fermions achieved by the 
time-dependent mass term (ma) 2 is typically smaller 
than the Pauli bound, nk = 1, even at the end of the 
modulus-driven phase. If one introduces the Yukawa cou- 
plings between two scalar fields <p, a and the fermion ip 
such as hi^tp^ 1 an d h2cnpip, the effective mass of fermions 
is expressed as [ETf 



TO e ff = TO + hl(f> + hiG. 



(4.1) 



In this case it is known that particle creation is most 
efficient when TO e ff vanishes, leading to rik ~ 1 both in 
the context of inflation Q and preheating [|23|j2j]-|27|] . In 
the present scenario, since m c Q vanishes or becomes close 
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